Study of pressure effect on the magnetic penetration depth in MgB 2 
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A study of the pressure effect on the magnetic penetration depth A in polycrystalline MgB2 was 
performed by measuring the temperature dependence of the magnetization under an applied pressure 
of 0.15 and 1.13 GPa. We found that A -2 at low temperature is only slightly affected by pressure 
= 1.5(9)%], in contrast to cuprate superconductors, where, in the same range of pressure, a 
very large effect on A~ 2 was found. Theoretical estimates indicate that most of the pressure effect 
on A -2 in MgB2 arises from the electron-phonon interaction. 
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Shortly after the discovery of the superconductivity in 
MgB2 at 39 K, several investigations of the pressure de- 
pendence of the superconducting critical temperature T c 
were carried outm&2£ Indeed, the magnitude and sign of 
dT c /dp may indicate a way to rise T c at ambient pressure, 
and moreover, help to understand the superconducting 
pairing mechanism. So far, all these studies show that T c 
decreases with increasing pressure, with a rate depend- 
ing on the method and the pressure medium used. The 
first hydrostatic measurement of T c (p) up to 0.7 GPa, re- 
veals that T c decreases reversibly under hydrostatic pres- 
sure at the rate dT c /dp = - 1.11(2) K/GPa^ The be- 
havior of T c (p) with pressure in MgB2 was attributed 
to the pressure induced lattice stiffening (increase of the 
phonon frequency) rather than to the decrease in the 
electronic density of states N(Ei?), that is only moder- 
ately affected by pressured Comparison with theoreti- 
cal calculations supported the view that MgB2 is a BCS 
superconductor with moderately strong electron-phonon 
coupling^ 

Apart from T c , an other relevant superconducting pa- 
rameter is the magnetic field penetration depth A. In 
fact the so-called superfiuid density A~ 2 is related to 
the Fermi velocity and to the density of charge carri- 
ers, and its temperature dependence gives information 
on the symmetry and on the magnitude of the supercon- 
ducting gap. A study of pressure effects on A~ 2 (0) can 
give important informations on how the electronic de- 
grees of freedom are affected by lattice modifications and 
on the nature of the electron-phonon coupling. Indeed, 
in cuprates high temperature superconductors (HTS), a 
huge pressure effect on A _2 (0) was found, in particu- 
lar in YBa 2 Cu 4 8 (Y124)£ Part of this effect was at- 
tributed to the strong renormalization of the Fermi ve- 
locity, and therefore to the effective mass, due to the non- 
adiabatic coupling of the electrons to the lattice. This 
result is in agreement with the substantial oxygen iso- 
tope effect found in cuprate o 9 ' 10 ! 11 ' 12 ! 13 and can be in- 
terpreted in the framework of non-adiabatic theory of 
superconductivity^ 



In this paper we report measurements of the mag- 
netic penetration depth under pressure in polycrystalline 
MgB2. The temperature dependence of A~ 2 (T) was ex- 
tracted from the Meissner fraction / measured in a low 
magnetic field. A small pressure (p) effect on A _2 (0) was 
found [ A £-y = 1.5(9)%], for pressure ranging from 0.15 
to 1.13 GPa. Theoretical calculations of the pressure ef- 
fect on A(0) confirm the smallness of this effect [ ^_ 2 ~ 
1.4%]. These results contrast with the huge effect (~ 40 
%) on A _2 (0) found in the non-adiabatic cuprate super- 
conductor Y124 in the same range of pressure^. 

The MgB2 powder sample was prepared by solid state 
reaction in flowing argon. As starting materials we used 
Mg flakes and amorphous boron (Alfa Aesar). A pellet 
with starting composition Mg x !B 2 was placed in a BN 
crucible and fired in a tube furnace under pure Ar gas. 
The sample was heated for one hour at 600° C, one hour 
at 800°C, and one hour at 900°C. 

The sample was first ground and then sieved in or- 
der to obtain a small grain size i?<10 /im, needed for 
the determination of A from the Meissner fraction mea- 
surements. The hydrostatic pressure was produced in 
a copper-beryllium piston cylinder clamp, especially de- 
signed for magnetization measurements under pressure 
(see Ref |l5|). The sample was put in a teflon cylin- 
der and the pressure cell was then filled with Fluorinert 
FC77 as pressure transmitting medium. The Meissner 
fraction was calculated from low field magnetization mea- 
surements (field cooling) performed with a commercial 
Superconducting Quantum Interference Device. In the 
approximation of spherical grains with average radius R, 
the Meissner fraction can be related to the magnetic pen- 
etration depth via the Shoenberg formulaii& 

a) 

Any change in x/xo due to pressure can be attributed 
mainly to a change of A(T), rather than to a change of 
R, which is practically pressure independent. The sam- 
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FIG. 1: Field cooled (0.5 mT) normalized magnetization of 
MgB2 as a function of temperature in the vicinity of T c for 
pz (open circles), pl (filled circles), and pn (open squares). 
The inset show the full normalized magnetization data for pl 
and ph as a function of the reduced temperature t = T/T c . 
Note that some of the data points were dropped for clarity. 



pie was measured at low pressure pl = 0.15 GPa and at 
the highest pressure available p B = 1.13 GPa. The pres- 
sure values were determined from a preliminary calibra- 
tion measurement, where the same sample was measured 
in the pressure cell together with a small piece of lead. 
The pressure was detected by measuring the T c shift of 
lead (T c (0kbar)=7.2K). We found dT c /dp = - 1.24(5) 
K/GPa, in good agreement with previous results^ The 
real zero-pressure measurement could not be performed 
in the pressure cell, since, in order to seal the cell, at 
least a small pressure had to be applied. 

In Fig. n the temperature dependence of the normal- 
ized magnetization of MgB2 at pl and pu is shown in 
the vicinity of T c , together with the zero pressure pz 
measurement performed on the same sample in a quartz 
tube. The T c 's were obtained from the intercept of the 
linear extrapolations (see FigQ: T c {pz) — 38.73(3) K, 
T c {p L ) = 38.66(3) K, T c {p H ) = 37.43(3) K. The magne- 
tization curves shift systematically with increasing pres- 
sure towards lower temperature. In the inset of Fig^ 
we show the normalized magnetization as a function of 
the reduced temperature t = T/T c for pl and ph- The 
identical temperature dependences for both pressure val- 
ues indicate the absence of stresses in the sample due to 
pressure. 

From the magnetization the Meissner fraction was ex- 
tracted. Because of the unknown mass of the sample 
in the pressure cell, the Meissner fraction for the lowest 
pressure (p£=1.5 kbar) was normalized at low temper- 
ature with the value obtained by measuring the same 
sample in a quartz tube at zero pressure pz- The val- 
ues of A -2 were then calculated using the Shoenberg for- 
mula (Eq. QJ). The pz and pl data were normalized 
at low temperature to the absolute value of A -2 mea- 
sured with the muon spin rotation technique (/^SR) on 
the same sample at zero pressured 



We define the shift of A 2 between two different pres- 
sures pl (lower) and p B (higher) at a temperature T as: 



AA- 2 _ \-\ PH ,T)-\-\p L ,T) 
A- 2 A- 2 ( Pi ,T) 



(2) 



The temperature dependences of A -2 for pl and pu, 
are shown in Fig. [21 together with the pz data measured 
in the quartz tube. The A~ 2 (T) curves for pz and pl 
overlap perfectly for temperatures not too close to T c , 
where a small pressure shift is present. Unfortunately, 
we were not able to obtain reliable data below 7K be- 
cause of a large background signal due to lead impurities 
in the pressure cell. Therefore, while for the pz measure- 
ment in the quartz tube we show the data down to 2 K, 
for the pl and the p B measurements we show the data 
only down to 7 K (see lower panel of Fig. 0) . Due to the 
pressure effect on T c , the curves for pl and pn are clearly 
distinct at high temperature, but they merge at low tem- 
perature. The error bars on A -2 (not shown) were deter- 
mined considering that magnetization measurements are 
reproducible within « 0.5 %. 

The shift of A~ 2 at the lowest temperature available 
T ~ 7 K between p L = 0.15 GPa and p H = 1.13 GPa 
(see Eq. EJ is: AA _2 /A~ 2 = 1.2(1.1)%. This result shows 
that only a very small pressure effect on A~ 2 is present. 
The pressure results are summarized in Table [I] 

To determine the pressure effect on the zero tempera- 
ture A _2 (0), a fit to the experimental data is needed. It 
was well established that MgB2 is a two band supercon- 
ductor with two superconducting gaps of different size, 
the larger one originating from the 2D c-band and the 
smaller one from the 3D 7r-bandi£ii2i2Si2i. Taking this 
into account, the temperature dependence of A, in the 
low temperature range, can be written in the foniS: 



A" 2 (T) = A" 2 (0) 



(1 - w) 



1 — w 



/ 2ttAi(0) 
V k B T 



1/2 



exp 



Ai(0) 
k B T 



2^A 2 (0) 
k B T 



1/2 



exp 



A 2 (0) 
k B T 



(3) 



Here, Ai and A 2 are the zero temperature small and 
large gap associated to the tt and a bands, respectively, 
and w = A^ 2 (0)/A- 2 (0). 

The zero pressure data were then fitted to Eq. |3]up to 
22 K. The fit gave the following results: A" 2 (0) = 91.1(2) 
fj, m- 2 , w = 0.22(3); Ai = 2.3(2) meV, A 2 = 7.4(4) meV. 
The estimated gaps arc in good agreement with previous 
results from penetration depth measurements 2 ^^. 

Because of the lack of data below 7 K for p^ and p B , 
these data were analysed with fixed w obtained from the 
fit to the zero-pressure data. Here it is assumed that 
w is not affected by pressure up to 1.13 GPa. This as- 
sumption will be shown below to be correct by our model 
calculations [see Eq.©]. The other fitting parameters 
(A _2 (0) and the gaps) were left free. As shown by the 
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FIG. 2: a) A~ 2 as a function of temperature for pz — GPa 
(open circles), pl = 0.15 GPa (filled circles), and pn = 1.13 
GPa (open squares), b) Low temperature region on a larger 
scale. The error bars (not shown) for temperatures below 10 
K, are about two times the symbol size and decrease with 
increasing temperature. The solid and dashed lines are fits to 
the pl and pn data, respectively, using Eq. [3] Note that the 
data points are much more dense than what is shown in this 
figure, since part of them were dropped for clarity. 



solid and dashed lines in panel b) of Fig. [2J the data 
are well described by Eq. ©. In particular the fitted 
curve of the p L data (solid line) follows very well the pz 
data below 7 K, as expected. The fit yields for A~ 2 (0): 
A- 2 (0) PL = 90.9(5)/otx" 2 and A~ 2 (0) Pff = 92.3(6)mto~ 2 
(see Table HP. So that the relative shift of A -2 at T = 
between p L = 1.5 kbar and pn = 11.3 kbar is [see Eq. 



AA- 2 
A" 2 



= 1.5(9)% 



(4) 



This result indicates that there is a small, but not zero, 
pressure effect on the magnetic penetration depth at 
K in the range p = 0.15 - 1.13 GPa. A summary of the 
pressure results are reported in Table |I] All the fitting 
parameters are summarized in the TablcllTI The two gaps 
show a small decrease with increasing pressure, in agree- 



TABLE I: Summary of the experimental and theoretical esti- 
mations of the pressure shift of A -2 . 





T 


PL 


Ph 


A" 2 (Pff ,T)-\-' 2 (p L ,T) 




(K) 


(GPa) 


(GPa) 


% 


exper. 


~7 


0.15 


1.13 


1.2(1.1) 


fit 





0.15 


1.13 


1.5(9) 


theory 





0.15 


1.13 


~ 1.4 t- 1.7 



ment with the corresponding decrease of T c . In Table ITT1 
we added also the ratio 2A/ksT c calculated for Ai and 
A2 at each pressure. This ratio results to be pressure 
independent within errors. 

Considering that pn — Pl — 1 GPa, Eq. J3J indi- 
cates that din X~ 2 (0)/dp is of the order of few %/GPa. 
Since A -2 (0) is proportional to the squared plasma fre- 
quency, uj 2 , then a free electron gas estimate would give 
din \- 2 (0)/dp = 1/B ~ 0.65%/GPa, where we have used 
lo 2 p oc I/O, being B = -dp/dhxQ, ~ 155 GPa the bulk 
modulus of MgB2f& and O the volume of the unit cell. As 
we show below, an improved estimate which takes into ac- 
count MgB2 band structure effects does not change much 
the free electron gas result. Let us consider the zero tem- 
perature expression of the penetration depths 



A" 2 (0) 



3ft7T 2 C 2 ^ 



ds\v n (s)\, 



(5) 



where the integral runs over the Fermi surface Sp of 
the n-th band (n = a, 7r) and v„(s) is the correspond- 
ing surface bare electron velocity vector. Let us use the 
model of Ref. |2(] where the it bands are modeled by a 
half-torus Fermi surface of area Si and Fermi velocity 
vi, while the a bands are approximated by a cylindrical 
Fermi surface of area Si and Fermi velocity V2 ■ Equation 
© then reduces simply to A _2 (0) = A^ 2 (0) + A^O), 
where A7 2/ 



'(0) oc SiVi, i = 1,2, and 
dlnA- 2 (0) 



dp 



where w = X 



dlnAr 2 (0) , .dlnAo" 2 ^) 
w J— — + (1 - w)- 



dp 

- 2 (0)/A- 2 (0) 



dp 



(G) 



SiVi/(SiVi + S2V2) is the 
parameter introduced in Eq.(|3J). Within the same ap- 
proximations, the electron density of states at the Fermi 
level reads Np = Ni + N%, where Ni oc ttSi/vi is the 
partial density of states for the band i = 1, 2. Therefore, 
since A~ 2 (0) oc S^i oc QS 2 /Ni, and considering that Si 
and S2 scale as Q~ 2 ^ 3 , Eq.© reduces to: 



dlnA" 2 (0)_ 1 
dp _ 3B 



1 — w + ryw 



T] + (1 - rj)N 2 /N F 



dhiNp 
dp 



(7) 



where we have introduced the parameter 77 = 
(d\nNi/dp)/(d\nN 2 /dp) whose calculated value ranges 
between r\ ~ Ij21 and r\ ~ Oi 2 ^ Hence, by setting 
d\nNp/dp ~ -0.31%/GPa^ to = 0.22, and N 2 /N F ~ 
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TABLE II: Fitting parameters obtained from the fit of the experimental data showed in Fig. |5]by the Eq. 



p 


A- 2 (0) 


Ai(0) 


Aa(0) 


w 


2Ai(0)/fc B T c 


2A 2 (0) /k B T c 


GPa 


(/im- 2 ) 


(meV) 


(meV) 











91.1(2) 


2.3(2) 


7.4(4) 


0.22(3) 


1.38(12) 


4.43(24) 


0.15 


90.9(5) 


2.4(1) 


7.3(1) 




1.44(6) 


4.38(6) 


1.13 


92.3(6) 


2.3(1) 


7.0(1) 




1.43(6) 


4.34(6) 



0.4,22 from Eq.Q we obtain din \- 2 (0)/dp ~ 0.8%/GPa 
or 0.5%/GPa according to whether rj = or r\ = 1, re- 
spectively. By the same token, it is easy to show that: 



dw 



w(l - w)(l - 77) d\nN F 



dp 7] + (1 - r))N 2 /N F dp 



(0-r -0.1)%/GPa, 
. . (8) 

confirming the assumption dw/dp = used in obtaining 
Eq.©. 

Our estimate AA^/A" 2 ~ (0.5-=- 0.8)% does not devi- 
ate from the free eletron gas result, suggesting that other 
factors than band structure should be considered in or- 
der to explain the measured value reported in Eq.©. 
It is then natural to consider the electron-phonon inter- 
action \ei-ph- The pressure effect on A e /- p /i in MgB2 is 
mainly due to a pressure induced hardening of the optical 
phonon modes. The electron-phonon renormalized pen- 
etration depth is A* _2 (0) = A~ 2 (0)/(1 + X e l- P h), where 
A _2 (0) is the bare quantity we have considered before. 
Hence: 



rflnA*" 2 (0) 
dp 



(0.5-9-0.8)%/GPa- 



A P _ 



ph 



d In A P 



-ph 



1 + K-ph dp 

(9) 

-1.7%/GPa^ and 
i (1.4^ 1.7)%, 

in better agreement with Eq.©. This simple analysis 
evidences therefore that the electron-phonon interaction 



which, by using dlnA e _ p h/ dp ~ 
Xel-ph — lr^ corresponds to AA~ 2 /A 



provides the main contribution to the pressure effect on 
the zero temperature penetration depth. This theoreti- 
cal estimate of the pressure shift of A(0) -2 is reported in 
the Tablc[I]for a direct comparison with the experimental 
finding. 

Similar experiment on the Y124, in the same range of 
pressure, gives a very large effect on A(0)~ 2 (~ 40%)£ 
Arguments are used there to deduce that part (~ 30%) 
of this large effect is due to the pressure dependence of 
the Fermi velocity, because of the non-adiabatic electron- 
lattice coupling. From the results of the present work, 
we can therefore argue that such non-adiabatic effects 
on A(0) -2 are negligible in MgP>2, as previously demon- 
strated by /zSR measurements^ 

In summary, we studied the pressure effect on the 
magnetic penetration depth at low temperature in poly- 
crystalline MgP>2. We found that pressure up to 1.13 
GPa induces a small positive change in A(0) -2 [ ^-2 = 
1.5(9)%], which is suggested to be due mostly to a pres- 
sure change of the electron-phonon coupling. 
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